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Abstract 

Within the framework of a model representing the potential 
of a graphene sheet U(z) as an electro-neutral layer formed 
by smeared carbon atoms, the effect of this potential on 
spectral characteristics of atoms adsorbed on a graphene 
sheet has been studied. Since the distance between the 
adsorbed atom nucleus and sheet surface significantly 
exceeds the radii of inner atomic shells, the potential U(z) 
makes influence on the continuum wave functions only. 
Their behavior in the upper semi-space (z>0) and in the 
lower one (z<0) where the adsorbed atom is located is 
defined by a jump of the logarithmic derivative of the wave 
function for z=0. The photoelectron angular distributions 
have been calculated for different mutual positions of the 
polarization vector e and the axis Z normal to the sheet 
surface. It has been shown that the existence of the electron 
waves reflected from the potential U(z) leads to evident 
asymmetry of the angular distribution relative to the plane 
2=0. The experimental observation of this effect is of great 
interest for photoelectron spectroscopy of atoms localized on 
graphene structures. 
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Introduction 

The number of papers devoted to investigation of 
carbon nanostructures and their feasible applications 
has seen rapid increases [see Lehtinen et al. (2003), 
Saito et al. (1992), Ajayan (1999), Gerber et al. (2010), 
Xia et al. (2008), Mishra et al. (2012), Rangel & 
Seminario (2010) and references therein]. The atomic- 
scale understanding of the properties of the adatoms 
on a graphene sheet has become an essential part of 
such studies. The photoelectron spectroscopy is an 
effective precise method to investigate the electron 
structure of adatoms [Carlson (1975)]. The study of the 
interaction processes of ionizing radiation with single 
adatoms and surface structures on their basis will not 
only provide an answer to a number of the intriguing 


fundamental scientific problems of atomic, molecular 
physics and surface physics but also promote the 
development of surface-structures nano-engineering. 

In this paper, the angular distribution of 
photoelectrons knocked-out from the deep subshells 
of atoms adsorbed on a graphene sheet was studied. 
Since we deal with physical adsorption when the 
adatom electronic structure is not subject to significant 
changes, the electronic subsystems of graphene sheet 
and adatom can be considered independently. The 
distance d between the adsorbed atom nucleus and 
sheet surface is several atomic units (1.87A [Lehtinen 
et al. (2003)]) and so it significantly exceeds the radii of 
inner atomic shells. That allows considering the wave 
functions of these shells with great accuracy as 
coinciding with the wave functions of the free atom. 
As far as the continuum wave functions are concerned, 
for near-threshold photoionization the wavelength of 
electrons produced in this process exceeds by far the 
bond lengths between the carbon atoms in the 
graphene sheet. This effectively obliterates the 
influence of the individual carbon atoms on the 
knocked-out low-energy photoelectron, thereby 
making the phenomenological potential U(z) formed 
by the z=0-plane-smeared carbon atoms to be an 
excellent approximation. A matrix element of the 
optical transition is formed near the adatom nucleus 
where the potential of graphene sheet is equal to zero. 
Therefore the specific shape of the function U( z) has 
no influence on the photoionization amplitude. This 
allows the real interaction potential U(z) having a 
finite thickness to be replaced by a zero-range 
potential U (z) = - AS(z ) providing a necessary jump of 
the logarithmic derivative of the wave function at z=0. 

A similar method based on the so-called Dirac-bubble- 
potential U(r) = -A8(r - R) [Lohr et al. (1992)] for the 
C6ofullerene shell was developed in [Baltenkov (1999), 
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Dolmatov et al. (2004), Baltenkov et al. (2010)] where 
this potential was used to describe the photoionization 
of endohedral atoms M@C6o as well. As a matter of fact, 
the M atom photoionization on a graphene sheet is the 
photoionization of the same atom M in the system 
M@C 60 where the fullerene shell is expanded on a 
plane. The difference between the structures of 
graphene formed by 6-rings carbon atoms and the 
fullerene shell consisting of 5- and 6-rings plays no 
principal role. 

Introducing the 8 - potential in the wave equation, per 
se, is a suitable mathematical way allowing 
replacement of the real interaction potential having a 
finite radius of action by a zero-range potential 
providing the necessary boundary conditions for the 
wave function of a particle moving in this potential. 
This idea originates from the classical paper on 
deuteron photodetachment [Bethe et al. (1935)]. This 
way allows avoiding arbitrariness for selecting the 
parameters of the potential well U(z). This 
arbitrariness is restricted by a single parameter A 
defined in experiment. It is evident that, as any model 
description, the 8 -potential model has a limited field 
of application. Within this model, it is impossible to 
consider the processes for which the electronic 
structure of the carbon nanostructures plays an 
essential role. Such processes are, in particular, those 
of C6o molecule photoionization. For them critically 
important are the geometrical properties of the 
fullerene shell, such as its diameter and thickness 
[Madjet et al. (2008), Riidel et al. (2002)]. The fullerene 
cage cannot be also considered as a source of the static 
potential for description of photoionization of atom M 
endohedrally confined in molecule C6o. For correct 
description of this process for some photon energies, it 
is necessary to use the more sophisticated approaches 
for the fullerene electronic hull [McCune et al. (2009)]. 

The sequences taking into consideration of the atom 
photoionization processes the carbon nano-structures 
potentials U(r) = -AS(r-R ) and U(z) - -A8(z) are 
principally different, which is understandable already 
at the qualitative level. In the first case, the spherically 
symmetric potential leads to radical changes in the 
photoionization total cross sections where the so-called 
confinement resonances discovered in the recent 
experiments [Kilcoyne et al. (2010)] appear. The 
amplitude of these resonances relate to the number of 
backscattered waves emitted by the nearby C atoms 
equidistantly located from the M atom [Baltenkov et al. 
(2010)]. At the same time, the photoelectron angular 
distributions of endohedral atoms M@C6o remain 
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unchanged, i.e. they are described by the same linear 
combinations of trigonometric functions as in the case 
of the free atom M; and there are only changes in the 
coefficients at these function. In the second case 
because of non-sphericity of the potential U(z), the 
principal changes should be expected, first of all, in 
the photoelectron angular distributions while the total 
cross sections will be slightly distorted owing to the 
small (as compared to C6o) coordination number, i.e. 
the number of carbon atoms equidistantly located 
from the adatom nucleus. However, the question 
whether the changes in the shape of the angular 
distributions are observable and what their character 
is can be resolved by calculating the differential cross 
sections. In this paper, an answer to these questions 
can be achieved by considering the photoionization of 
atom near the graphene sheet within the framework of 
a simple model. 

Differential Photoionization Cross Section 

In the spherical coordinate system connected with the 
sheet surface, the axis Z is considered to be directed 
along the normal to the surface. The ionized atom M is 
placed in the lower semi-space (Fig. 1). 



FIG. 1 THE SCHEME OF THE ADSORBED M ATOM LOCATION 
ON THE GRAPHENE SHEET. THE PLANE WHERE THE ATOM 
CARBONS ARE LOCATED IS PERPENDICULAR TO THE PAGE 

The spherical coordinates of the radiation polarization 
vector e are (9 e , </> e ) and the coordinates of the vector 
of photoelectron momentum k are ( k , 9, <j>) . It is 
known that the differential cross section of dipole 
photoionization of atom for linearly polarized light is 
determined by the cosine of the angle 9 ek (between 
the vectors e and k) and has the general form (Bethe, 
H. 1933) 

d °- = ^-[\ + pP 2 (cos3 ek )] (1) 
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Here a is the total photoionization cross section; // is 
the dipole asymmetry parameter. In the above- 
described coordinate system, the argument of the 
Legendre polynomial P 2 (cos 9 ek ) is represented as 

cos9 ek =cos.9 e cosi9 + sin 1 9 1? sin l 9cos(^ e -</>) (2) 

For the fixed coordinates of the radiation polarization 
vector 9 e and <j > e , the shape of the electron angular 
distribution in the coordinate system connected with 
the sheet surface will be defined as the following 
function of spherical angles 9 and (f> 

> 1 . (3) 

Since we are interested in the shape of the angular 
distribution only, this function was normalized so that 
it is equal to unity for 9 ek = 0 . The vector J = k F(9,</>) 
is a flux of photoelectrons knocked-out by the ionizing 
atom (here and throughout the text the atomic system 
of units e=m = h = l is used). The Cartesian 
components of this flux are defined by the following 
formulas 


j x = k x F (9, <f>) = kF (9, </>) sin 9 cos </> , 
j y =k y F{9,</)) = kF{9,<l))s\n9sin<l) , (4) 

j. = k z F(9, <l) = kF(9, <!>) cos 9 . 

The photoelectron wave function that corresponds to 
the flux J was represented as 

¥k (r) = F(9,<P)e lk r . (5) 

The electron wave function in the form (5) describes 
the photoelectron as a free particle since the Coulomb 
field of the atomic residue gives for the electron flux 
density the corrections that go to zero at — > oo 
[Landau & Lifshitz (1965)]. 

Graphene Sheet Effect 

Until now, the above-given formulas relate to the free 
M atom. In the case of the ionization of adsorbed atom, 
the potential plane formed by the graphene sheet 
distorts the final-state wave function (5). Near the 
photoeffect threshold when the photoelectron wave 
length A significantly exceeds the distance between 
the carbon atoms in the 6-rings forming the graphene 
sheet, its potential (by analogy with the Dirac-bubble- 
potential for C6o shell [Lohr et al. (1992)] was regarded 
as a phenomenological potential U(z) = -AS(z) to 
simulate the smeared carbon atoms. 

Within this model description, the wave function of 
photoelectron leaving the atomic residue will obey the 
following Schrodinger equation 


A y/ k +2 [E-U{z)]y/ k =0; 

E = k 2 1 2 = (k 2 + k 2 + k 2 ) / 2 . (6) 

The solution of this equation is written as 

¥k ( r ) = F(94)e Kx+lKy 4>(z). (7) 

It is evident that for U(z) = 0 the function (7) coincides, 
as it should be, with the function (5). For the function 
t/>(z) we have the equation 


~Y + [^z + 2A5(z)]</>(z) = 0 . (8) 

oz 

The solution of Eq.(8) we write in the form (Flugge, S. 
1971) 

(j){z) = e lk ‘ z + Be~ ,k ’ z for z<0, 

<p(z) = (1 + D)e ik ’ z for z>0. (9) 

In the formulas (9) three waves are distinguished: an 
incident wave with unity intensity, a reflected wave 
with intensity j B | 2 and a wave passed over the 
potential well with intensity 1 1 + D | 2 . Because of 
continuity of the function </>(z) at z=0, the relation B=D 
for the scattering amplitudes was obtained. 
Integrating the equation (8) over z within the limits 
±s — » 0 , we obtain the equation connecting a jump of 
the logarithmic derivative of the wave function 
A L(z = 0) with the strength of the delta-potential: 


AL = 


dt/t , 

dtp , 

-N I +£■ 

. dz 



1 

m 


= -2 A . 


( 10 ) 


Substituting the functions (9) in Eq.(10), we obtain the 
following expression for the scattering amplitudes 


B = D = - 


A 

A + ik z 


( 11 ) 


Owing to electron wave reflection from the potential 
plane, the electron flux vector has the following z- 
component in the lower semi-space 


j[=F{9,m+\BUk^s9, ( 12 ) 

and in the upper one 

j“ =F(9,</>)\l + D\ 2 kcos3 = [l-\B\ 2 ]kcos9. (13) 

As far as the Cartesian components of the 
photoelectron flux j x and j v are concerned, they are 

undistorted by the potential U(z) and defined as 
before by the formulas (4). Finally, the formulas 
associating the function F(9,<l > ) for the free atom with 
the angular distribution of electrons knocked out of 
the atom localized near the potential well U (z) in the 
Cartesian coordinates have the following form 


S x = F(9,<f>)sin9cos<f> , 
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S v = F(9,<f>)sin9s'm<f> , (14) 

S. = F(9,<f>)[\- 1 B | 2 ]cosi 9, for the upper semi-space, 

S. = F(9,<f>)[ 1+ 1 B | 2 ]cos6>, for the lower semi-space. 
The function \B\ 2 in (14) has the form 


\B\ 2 =- 


2 — 2 2 ~- ( 15 ) 

A 2 +k~ cos 2 9 

From (15) it follows that when the electron moves 
from the upper semi-space to the lower one or vice 
versa, i.e. at angles 9 = n 1 2 or 3^/2, the z- 
components of the function S z Eq.(14) have the jump 
due to partial reflection of the electron flux from the 
potential well. 


Let us investigate how the derived formulas are 
modified if the finite thickness of the graphene film is 
taken into consideration. The formulas (9) will have 
now the following form 

<f>(z) = e lk * z + Be~ ,k * z for z<-b/2, 

<t>(z) = A^(z) + C£(z) for b/2>z>-b/2, (16) 

</>(z) = (\ + D)e ik ~- z for z> b/2. 

Here A%(z) + CC(z) is the linear combination of the 
regular and irregular solutions of the Schrodinger 
equation inside the graphene sheet. Matching the 
logarithmic derivatives at the points z = ±b / 2 (here b 
is the effective film thickness shown in Fig. 1 by the 
dashed lines), the following expression for the 
amplitude of the reflected wave was obtained 




L_+L + 


(17) 


Here L + and L are the logarithmic derivatives at the 
points z = ±b / 2 , i.e. the two model parameters to be 
defined. Let us suppose that the potential U(z) is a 
squared well with the depth W and the thickness b. 
Then the linear combination of the solutions within 
the range b/2>z>-b/2 can be represented as 


A%(z) + C£(z) = A sin Kz + Ceos Kz , (18) 

where the wave vectors K and k 7 are connected with 
the relation K 2 / 2 = k 2 12 + W . The amplitude of 
electron wave reflection within this model is defined 
by the following formula 

i^i 2 = — — ~ k ?\ — 2 — • ( w ) 

1 1 (k 2 +K 2 ) 2 + 4k 2 K 2 cot 2 Kb 

The amplitude of the electron wave passing over the 
potential well 1 + £>| 2 is connected with (19) as follows 


l-\B\ 2 = |l + Z)| 2 . (20) 

As before, the formulas (14) describe the shape of the 
photoelectron angular distribution and they will be 
used for numerical calculations. 

Numerical Calculations 

S-atomic State 

For the strength of the delta-potential A defining the 
scattering amplitude (11) to be evaluated, it was 
assumed that the potential of the graphene sheet U (z) 
is similar to a plane-expanded potential of the 
fullerene C6o shell. Based on this assumption, the 
parameter A can be found from the experimental data 
on the fullerene radius R and the electron affinity 
energy I. According to [Baltenkov (1999)], 

A = j3(l + cothj3R)/2, (21) 

where the parameter /? is /? = (2 / ) 12 . For I = 2.65 eV 
and R = 6.639 [Tossati et al. (1994)] the value of the 
potential strength A is 4=0.441 . 

This value corresponds to spherical potential well of 
the squared profile with the depth W and the 
thickness b that are connected with each other so that a 
bound state with energy I exists in this potential well. 
In the case of C6o, it is generally accepted that the 
fullerene electronic hull has the thickness b=\ .89 and 
the depth W =8.24 eV [Baltenkov et al. (2010)]. It is 
evident that for the graphene sheet these values of the 
parameters A, W and b should be considered as 
tentative. At the end of this section, we will analyze 
how the shape of photoelectron angular distribution 
are modified with variation of these parameters. 

Let us consider the case of the s-states of the ionized 
atom. For them, the dipole asymmetry parameter is 
P = 2 and the function F(9,(/)) = cos 2 9 ek . The 
coordinate system is chosen so that the X-axis is in the 
same plane as the vectors e and Z where the intensity 
of photoelectron emission is maximal. In this 
coordinate system, the azimuth angle (f> in the 
formulas (14) is equal to zero and the photoelectron 
spectrum shape is defined by the x- and z-components 
of the function S. The calculation results for the 
function S defining the photoelectron spectrum shape 
for two polar angles of the e vector 9 e = n / 2 and 0 are 
presented in Figs. 2 and 3. The angular distributions 
for free atom photoionization are presented in these 
figures and in the other ones for comparison. 

In Fig. 2 the polar angle 9 e is 9 e = n / 2 and the 
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function F(9,</>) is F(9,</>) = sin 1 9 . This function is not 
equal to zero at the angles 9 = n / 2 and In 12 and at 
these points the jumps of the function Si manifest 
themselves in the explicit form. The upper spectrum 
part corresponding to the angular distribution of 
photoelectrons passing through the potential plane is 
evidently less than the lower one. With the rise in 
photoelectron kinetic energy s the reflected wave 
amplitude (15) decreases and the spectrum shape for 
adatoms is transformed into that corresponding to the 
free atom. 

The curves for the angle 9 e = 0 for which the 
function F(9,</>) is F(9,</>) = cos 2 9 are presented in Fig. 
3. For the angles 9 = n/2 and 3zr / 2 the 
function F(9,t/>) vanishes and so the jumps in Sz are 
not explicit. The upper spectrum part corresponding 
to the angular distribution of photoelectrons passing 
through the potential plane, as in Fig. 2, is vividly less 
than the lower one. For e || Z the spectrum shapes are 
independent of the azimuth angle <j) and the 30- 
pictures of the photoelectron angular distribution are 
the figures formed by rotating the curves around the 
Z-axis. For all the polar angles of the vector e, the 
picture of spectrum transformation due to the 
potential U (z) is qualitatively the same: electron wave 
reflection leads to enhancing the flux to the lower 
semi-space and weakening to the upper one. With the 
rise in photoelectron kinetic energy e the adatom 
spectra are transformed into those corresponding to 
the free atom. 

1.0 


0.5 


0.0 


-0.5 


- 1.0 
-1 

FIG. 2. THE ANGULAR DISTRIBUTIONS OF PHOTOELECTRONS 
KNOCKED OUT OF THE S-ATOMIC SUBSHELLS FOR 
DIFFERENT ELECTRON MOMENTUMS K. THE POLARIZATION 
VECTOR IS e || X . DASH-DOT-DOT BLACK LINE - FREE ATOM; 

SOLID RED LINE - ATOM + GRAPHENE SHEET FOR K=0.6; 
DASHED BLUE LINE - ADATOM FOR K=1 .0; DOT GREEN LINE - 
ADATOM FOR K= 2.0. 
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FIG. 3. THE ANGULAR DISTRIBUTIONS OF PHOTOELECTRONS 
KNOCKED OUT OF THE S-ATOMIC SUBSHELLS FOR 
DIFFERENT ELECTRON MOMENTUMS K. THE POLARI Z ATION 
VECTOR IS e || Z . DASH-DOT-DOT BLACK LINE - FREE ATOM; 

SOLID RED LINE - ATOM + GRAPHENE SHEET FOR K= 0.6; 
DASHED BLUE LINE - ADATOM FOR K=1.0; DOT GREEN LINE - 
ADATOM FOR K= 2.0. 

Let us compare the amplitudes of electron wave 
reflection in the models of zero and finite thicknesses 
of the potential well U (z) . In the 8 -potential model 
the reflection coefficient is defined by the formula (15). 
The value of the potential strength A in (15) 
corresponds to potential well with the above defined 
parameters W and b. Another several pairs of the 
parameters W and b meeting the condition of bound 
state existence in the well [Baltenkov et al. (2010)] were 
taken into consideration. The comparison results of the 
reflection amplitudes (15) and (19) with these 
parameters are presented in Fig. 4. It is seen that the 
narrower and deeper the potential well is, the wider 
the electron energy range where the behavior of the 
reflection amplitudes in the both models is similar. 



FIG. 4. THE AMPLITUDE OF ELECTRON WAVE REFLECTION IN 
THE 8 -POTENTIAL MODEL AND FOR SOME SQUARED 
POTENTIALS. 

In Fig. 5 it is possible to see the spectrum 
transformation with changes in the strength of the 
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delta-potential A for fixed electron energy (k= 0.6). 
With the rise in the parameter A, the reflected wave 
amplitude (15) goes to unity and the photoelectron 
flux is almost totally reflected by the potential plane to 
the lower semi-space. 



FIG. 5. THE ANGULAR DISTRIBUTIONS OF PHOTOELECTRONS 
KNOCKED OUT OF THE S-ATOMIC SUBSHELLS FOR 
DIFFERENT VALUES OF THE PARAMETER A AND FOR FIXED 
ELECTRON MOMENTUMS K= 0.6. THE POLARIZATION VECTOR 
IS e || Z . DASH-DOT-DOT BLACK LINE - FREE ATOM; SOLID 
RED LINE - ATOM + GRAPHENE SHEET FOR A=0.441; DASHED 
BLUE LINE - ADATOM FOR A=1.0; DOT GREEN LINE - 
ADATOM FOR A=2.0. 

Non-zero Orbital Moments 

The photoelectron angular distribution for the atomic 



shells with nonzero orbital moments is characterized 
by the dipole asymmetry parameter /? changeable 
within the range -1 < /? < 2 . It is possible to see the 
evolution of the photoelectron spectrum shape for the 
parameter p in the interval -0.5 < // < 1 .0 in Fig. 6. 
The polarization vector e here is parallel to the Z-axis. 
The calculations were performed for two values of 
electron momentum k= 0.6 and k= 1.0. As before, with 
the rise in electron energy e the adatom spectra are 
transformed into those corresponding to the free atom. 
The electron angular distribution for parameter p = 0 
corresponds, per se, to the picture of the light intensity 
of the isotropic source (candle) located near 
semitransparent mirror when an observer is far from it 
at the distance significantly exceeding d and so the 
light source and its reflection in mirror are perceived 
as a point source of light. 

As it follows from the figures, the existence of the 
electron waves reflected from the graphene sheet 
potential U(z) leads to quite observable asymmetry in 
the angular distribution relative to the plane z =0; the 
less the photoelectron kinetic energy s is, the more 
pronounced the effect is. 




s. 
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Adatom k=0.6 

Adatom k=1 .0 

s x . 


FIG. 6. THE ANGULAR DISTRIBUTIONS OF PHOTOELECTRONS KNOCKED OUT OF THE ATOMIC SUBSHELLS WITH NONZERO 
ORBITAL MOMENTS FOR DIFFERENT DIPOLE ASYMMETRY PARAMETERS P AND DIFFERENT ELECTRON KINETIC ENERGY 

£ = k 2 / 2 . THE POLARIZATION VECTOR IS e 1 1 Z . DASH-DOT BLACK LINES - FREE ATOM; SOLID RED LINES - ATOM + 
GRAPHENE SHEET FOR K= 0.6; DASH BLUE LINES - ATOM + GRAPHENE SHEET FOR K=1.0. THE 3D-PICTURES OF THE 
PHOTOELECTRON ANGULAR DISTRIBUTION ARE THE FIGURES FORMED BY THE CURVES ROTATION AROUND THE Z-AXIS. 
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For k= 0.6 the electron wave length is A, = 2tt / k « 1 0 
while the distance between the carbon atoms in 
graphene is a=1.42A=2.68 [Abergel et al. (2009)]. For 
these electron wave lengths A» a the developed 
model is parametrically justified The experimental 
observation of this new physical effect is of great 
interest for photoelectron spectroscopy of atoms 
localized on graphene structures. 
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